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Abstract 

The peak effect lias been investigated in clean Nb crystals with artificially 
corrugated surfaces by measuring the linear surface impedance in the IkHz- 
IMHz frequency range. From a two-mode analysis of the complex spectra, 
we establish that vortex dynamics is governed by surface pinning and deduce 
the associated vortex slippage length. We demonstrate experimentally and 
theoretically that the peak effect is related to a transition from collective to 
individual surface pinning. A proper account of the peak-effect anomalies im- 
plies softening of the shear rigidity by disorder-induced lattice deformations. 

This leads to a vortex crystal-glass transition induced by surface defects. 
PACS numbers: 74.60.Ge 
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The peak effect (PE) is a well known anomaly of the mixed state due to vortex lattice 
(VL) softening in the presence of random defects [1,2]. It is observed at the upper critical 
field Bc2 of clean superconductors [3,4] , well below Bc2 in disordered oxide superconductors 
(second peak), or at the vortex-melting fine in very clean YBa2Cu307 crystals [5]. The 
physics of PE is the order- disorder transition, which is of interest for a broad range of elastic 
media. Recent developments concern vortices in conventional materials (Nb, NbSe2) , where 
PE is investigated in great detail [6-8] . Originally believed to be a crossover from collective 
to individual vortex pinning according to the Larkin-Ovchinnikov (LO) scenario [2], PE is 
now considered as a genuine transition [6,8], with a quasi-long-range order (Bragg-Glass) 
in the low-field phase [9,10]. However, strong limitations for a quantitative analysis of the 
phenomenon arise from the poor knowledge about the random defects responsible for PE. 

In this work we consider a new situation where this difficulty can be overcome : the 
peak effect due to the surface pinning of vortices. The pinning landscape is provided by 
the surface corrugation C(r) of a slab in perpendicular or tilted magnetic field; a controlled 
roughness Q* ~ 1-lOnm is produced by ion beam etching (IBE) and measured by atomic- 
force microscopy ( AFM) . Far from the upper critical field 3^2 the intervortex interaction is 
strong and suppresses the surface pinning. This means that the surface pinning is collective 
similar to the collective pinning in the bulk [2]. But close to Bc2 the VL distortion by 
surface pinning become essential, and near the surface the ordered VL transforms to the 
disordered vortex array, which we call surface vortex glass. Formation of the surface vortex 
glass is accompanied by the crossover from collective to individual pinning. This results in 
the growth of the critical current, i.e. in PE. 

Experiments are performed on a series of millimeter-thick slabs of very pure Nb. Starting 
from electro-chemically polished samples, a controlled corrugation is introduced by ion- 
beam etching (IBE) the surfaces with a 1.5mA/cm^ fiux of SOOeV-Ar"*" ions for 10-100 
min. Given the sputter yield, 0.6 Nb atom per Ar"*", and the atomic spacing c = 0.26 nm, 
the average etching rate is of the order of 100 nm/min. From this stochastic process, one 
expects a roughness amplitude C,* ~ (cA^;)^/^ and a white spectrum S'^(k) = / dr e~*''""(C(r-|- 
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R)C(R))r = (^IS.zj'K in the range — c"^ of wave number k. In practice, the wave number 
cut-off kc <'\-lc due crystal rearrangement associated with local heating; whence a reduction 
of by a factor /CcC^ 10~^. This picture is confirmed by AFM-inspection, which shows a flat 
spectrum below IjivcT^ and a power law dependence S(^{k) ~ l/k-l/k"^ above. The etched 
surfaces have a broad-band corrugation, in the range k — l-100/xm~^ covering the vortex 
reciprocal-lattice unit Q = Gq^ ~ (50nm)~^, which is very effective for VL-pinning. Here 



oo = ^Jipo/nB is the radius of the vortex Wigner-Seitz cell (on the order of the intervortex 
spacing) and (po = h/2e is the flux quantum. 

The pinning strength is deduced form the AC response. The AC surface impedance 
Z{u!) — —i^qujXac is determined by the effective penetration depth given by 

1 1 / 1 . ^ IsB 



Here ct/ is the flux-flow resistivity, s^Pq is the vortex- line tension (energy per unit length), 
and Ac the Campbell penetration depth for bulk pinning. The surface pinning is taken 
into account by the length L5 ~ 0.1-1 OO/^m, which can simulate the effect of bulk pinning 
characterized by the Campbell depth A^. This expression was derived within the frame of 
the two-mode electrodynamics [11,12], which incorporates the surface pinning by introducing 
a phenomenological boundary condition, 

/u(0) (9u(0)\ ^ 



imposed on the average VL displacement vl{z) = u(r, z) at the surface of the sample, which 
occupies the semispace 2; < 0. Here Is is a slippage length and the displacement u is averaged 
over the position vectors r in the xy plane. Physically this boundary condition presents the 
balance between the pinning force —etpQuiQi) /Is and the line tension force eipoduiO) / dz due 
to the vortex bending. 

Figure 1 shows the PE due to surface pinning. Data relate to sample #S2 with dimensions 
25 X 10.1 X 0.87 mm^. The sample was annealed in ultra-high vacuum which gives a low 
residual resistivity pn = H nQcm and 3^2 — 0.29 T at A.2K. This value is consistent with 



earlier measurements within the anisotropy of Nb crystals {Bc2{0) — 0.280-0.295 at 4.2K) 
[13]. Prom the 90 min IBE exposure we estimate C* ~ 5 nm, in qualitative agreement with 
AFM measurement. Data points are obtained by fitting the penetration-depth spectra as 
shown in the inset of the figure. Metastability in the vortex number and/or arrangement is 
removed by feeding a large transient current (7 ==40 A) in the sample prior to measurement. 
The abrupt onset of the AC-flux penetration along the samples edges which are parallel to 
the field precludes quantitative analysis for B ^ 0.95Bc2; this difficulty can be overcome 
by working in oblique field (figure 1). More experimental details including the calibrating 
procedure can be found in Refs. [12]. Already present in pristine samples, the PE is strongly 
enhanced by IBE up to a factor of 3 at the maximum dose. This contrasts with conventional 
mechanical or chemical treatments which show ineffective in increasing the PE in our Nb 
samples. This distintive property of IBE emphasizes the importance of small-scale atomic 
corrugation for the surface pinning of vortex lines. 

We separate quantitatively the bulk and surface contributions by relying on the full 
IkHz-lMHz spectrum XAcif) [12]- From fits of XAcif) with Eq.(l) we deduce Ls, Xc and 
(7/. Remarkably enough, we always find ^ 10m, meaning that the bulk pinning strength is 
neghgible. In fact there is no difference in the spectra taken on both sides of the peak which 
are equally well described by Eq.(l). This confirms that surface pinning is most relevant 
for the vortex dynamics in our experiment. The oblique- field data are larger by a factor 
~2; this general trend is due to surface- reinforcement of superconductivity in tilted fields. 
Similar data are observed at lower temperatures, with, however, larger Bpf; = 0.95Bc2 (1.8K) 
and a less pronounced PE. 

Using the Abrikosov expression, //o£— {Bc2 — B)/2.32k^ for B < Bc2, with «; = 1.3, and 
i?c2 = 0.29T, we deduce from Eq.(l) the ls{B)-data in Fig. 2. By taking the bulk expression 
for e, we disregard surface enhancement in oblique field and therefore overestimate the Ig^ 
data in this geometry. The high- field plateaus Is — h above Bpk in Fig. 2 are suggestive of 
an individual pinning, since in this case the slippage length is given by a typical curvature 
radius of the corrugated surface, which should not depend on the vortex density (magnetic 
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field). This interpretation is supported by the fact that the contact angle for VL at the 
surface ao//o — 0.1, estimated from Iq ~ 0.5/im (normal field), agrees with the maximum 
tilt C*/<^o^ 0.1 deduced from surface corrugation. This suggests that the surface-pinning 
strength above Bpk is indeed limited by geometrical considerations. By contrast, the strong 
suppression of 1/ls below Bp^ (factor ~ 10 in oblique field) refiects the collective regime of 
surface pinning, which was known earlier in rotating ^He [14] . The transition is sharp unlike 
the continuous ones reported in Refs. [6,8]. Thus the experiment provides an evidence that 
PE is accompanied by the crossover from the collective to the individual surface pinning. 

We start the theoretical analysis from calculation of the elastic response of the semi- 
infinite vortex lattice to the surface force presented by one Fourier component f(r) = 
f(k)e*'*'". Any surface force on a vortex should be balanced by the line-tension force: 
f(r) = S(podu/dz. Thus the force produces the displacement field in the bulk, which is 
assumed to be transverse (U(k) _L k), since the VL is much more rigid with respect to the 
longitudinal force than to the transverse one. The equation of the elasticity theory is: 

■C66A;2 + C44(k)A;,2]U(k) = 0, (3) 

where Cqq is the shear modulus and 

is the tilt-modulus, which takes into account nonlocal effects due to long-range vortex- vortex 
interaction. The solution of equation (3) yields two values for k"^ = — p^, which correspond 
to the two modes generated in the bulk. Usually Cqq <^ sB <^ B^/iiq and approximately 

2 (1 A^A;^)C66 , 2 2 I ;2-^ 2 /r\ 

By,o + eBX^k^ ^ ^^^^^^^^ 

The solution is a superposition of two components, U{x, z) — e^^^ {uic^^^ + U2e^^^), with the 
amplitudes ui and U2 determined from the boundary conditions. The first one is given by the 
surface force: /(k) = e^poiPiUi +P2U2), and the second one is that the tangential magnetic 
field determined by the London equation, h(k) — ikzB\J{k) [1 + A^(/c^ -|- kl)]~^, vanishes at 



the sample border. Eventually one can find the linear response relation /(k) = C(A;)C/o(k), 
which connects the force and the displacement U{x, 0) = U^e^^^ at the surface in the Fourier 
presentation, where Uq — ui + U2 and 



Most of the deformation energy is stored in the first long-wavelength mode. 

The linear-response function allows to calculate the random displacements produced by 
the surface pinning. Considering that vortex cores end perpendicular to the local surface, 
the corrugation produces a force /m(i'i) = —e(pod({ri)/dXjn applied to the end of the vortex 
line, where rj = Ri + Uq is the 2D position vector of the vortex end. For small displacements 
Uo, the balance of force yields the boundary condition 



aUo^(Ri) ^ aC(Ri) 



= . (7) 



dz dx^ 

In the Fourier presentation the random force is f(k) = — ik£(/?o / dr e~'^^C,{Y). Since C{k) oc k 
at small k, the integral for the mean-square-root displacement is divergent at small k. But 
we can obtain a convergent integral for the mean-square-root shear deformation 



;(VUo)') = {{dUoJdy + dUoy/dx) 

Here the integration over k is fulfilled over the Brillouin zone and Q is the recipro- 
cal vortex-lattice vector. We shall approximate the surface corrugation spectrum as 



S'(^(k) = 2n(*'^r'^e~''^'^, where (* = -y/(C(R)^)R is the amplitude of the surface roughness 
and Td is the correlation length of the random pinning force, which cannot be much less 
than the coherence length ^. In the collective regime k ^ Q, considering a broad corru- 
gation spectrum (r^ <^ oq), summation over Q can be replaced by integration, and we 
obtain 



where lo — r^/C* is on the order of the geometric curvature radius. Here we have used 
the large k hmit of Eq.(6), C{k) ^ kipo{CQQe / B)^/'^ , which is a good approximation when 
A > ao- 

In the experiment the AC fields produce additional quasistatic displacements u(Ri) 
superimposed on the static displacements induced by pinning. By differentiating Eq.(7) for 
small u, we obtain the boundary condition: 

In the Fourier presentation the uniform displacement (k = 0) is coupled with the nonuniform 
displacements (k 7^ 0). Treating the latter by the perturbation theory we arrive at the 
the boundary condition Eq. (2) imposed on the averaged, i.e., uniform displacement with 
slippage length given by 



Using the same approximation as in calculating ((VUo)^) we receive for the inverse slippage 
length 



Our derivation of the collective pinning strength oc I/Z5 was not based on the heuristic 
concept of Larkin-Ovchinnikov domain [2] usually used for the analysis of collective pinning. 
The domain size Lc is determined by the balance of the pinning energy and the elastic energy. 
The average energy density of pinning weakened by collective effects is Ep ~ e(pou'^ / alloy/Nc, 
where = Ll/a^ is the number of vortices in the Larkin-Ovchinnikov domain. In order 
to estimate the elastic energy one should know how far the shear deformation from surface 
disorder penetrates into the bulk. This length is not evident apriori, but aposteriori, from 
our analysis. Both experiment and theory show the importance of the short lengthscales, 
~ oo for the surface pinning; as a consequence, the deformation penetration length. 




X 



C{k) 



(11) 




(12) 
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~ ~ ysf^o/C'ee which corresponds to large k hmit in Eq. (5), is finite. Then the elastic 
energy per unit area is Eei ~ {Cqq/pi){u'^ / L^) ~ v^C'ee^^M^/L^. The condition Ep ~ E^i 
yields ~ IqUq^ Cqq/ SifQ and since Is ~ lo\/N'c — ioLc/ao one obtains for Is the value of 
the same order as in Eq. (12). 

Let us apply the received expressions for the magnetic fields far from Bc2- In this case 
e ~ {ipo/iioX'^)ln{Bc2/B), Cge ~ (fioB/iioX^ and according to Eq. (12) Is oc ^J B / ln{Bc2/ B). 
Then the surface length L, oc [B/\n{B^2/B)]^/'^ grows with the magnetic field in qualitative 
agreement with the experiment (Fig. 1). This is a regime of collective pinning when Is > lo- 
At the same time since <^ ao the vortex lattice shear deformation remains small according 
to Eq. (9). Close to Bc2 e ~ {Bc2 - B)/hqk^, Cee ~ {Bc2 — BY/k^, and ~ Oq ~ 
Then from Eqs. (9) and (12) we receive that ((VUo)^) ^ (^V^o)^c2/(5c2 - B) and Is ~ 
{IqI {Bc2 — B)/Bc2- Thus Is decreases when B approaches to Bc2 and becomes smaller 
than Iq for Bc2 — B < Bc2^'^/ll- This means that surface pinning ceases to be collective and 
the transition to the individual pinning occurs. At the same time, at Bc2 — -B ~ Bc2C'^/lo, 
the deformation ((VUo)^) becomes of the order of unity, which means destruction of the 
crystaUine order at the surface. The state without the crystalline order near the surface can 
be called surface glass. Thus the crossover from the collective to the individual pinning is 
accompanied by the crystal-glass transition at the surface. 

Now let us check that this transition, or crossover, can explain PE. Despite the fact that 
Is oc V Bc2 — B decreases at B approaching to Bc2, according to Eq. (l) Ls continues to 
grow proportionally to B^z — B whereas in the experiment (Fig. 1) Ls{B) decreases 
close to the peak. However, the growth of the deformation ((VUo)^), which accompanies 
the decrease of Is eventually invalidates the perturbative approach used above. Qualitatively 
this can be corrected by introducing the renormalized deformation-dependent shear modulus 
in the form Cqq — CQe{l — a ((VUo)^)) ~ Ceei^ — B/Bpk). Here the field Bp^ corresponds 
to the crystal-glass transition, where Cqq = 0, and a is an unknown numerical factor, which 
could be close to 0.1 as in the Lindemann criterion. Using renormalized modulus Cqq in place 
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of Cee in Eq. (12) we obtain that Is as well Ls [see Eq. (1)] decrease inversely proportionally 
to ^ Bpk — B in qualitative agreement with experiment (Figs. 1 and 2). On the right of the 
peak, where pinning is individual and the shear rigidity is absent, Is ~ h does not depend 
on the magnetic field B, while Ls oc l/{Bc2 — B) grows with B. 

The close relation between PE and vanishing of the shear modulus of VL was suggested 
in the early studies of PE [1,2]. The new feature in our scenario is that the shear vortex 
elasticity vanishes only in some surface layer. The width ~ 1/pi oc 1 / V Cge of this layer 
diverges if Cee and finally surface disorder can contaminate the bulk of the sample. But 
this should happen for fields B > Bp^, since the peak field Bp^ is determined by Cqq. Our 
scenario agrees with recent STM imaging of the vortex array by Troyanovski et al. [15]. This 
experiment has revealed the disorder onset on the surface of 2H-NbSe2 sample, coinciding 
with the onset of the peak effect. Troyanovski et al. related the disorder onset with the bulk 
pinning. In order to discriminate two scenarios it would be useful to supplement the STM 
probing of the vortex array at the surface by probing vortex arrangements in the bulk. 

In conclusion, this experiment and the accompanying theoretical model provide a new 
plausible scenario for the peak effect. This is a comprehensive demonstration of a new 
situation: the order-disorder transition in a vortex lattice triggered by a weak surface disor- 
der. Beside its experimental relevance, this new established mechanism offers an interesting 
paradigm for elastic systems at the upper critical dimension for disorder. 
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FIG. 1. Peak-effect in tfie elastic response Lg (B) of a surface-pinned vortex array. Diamonds 
and full circles correspond to oblique and perpendicular field orientations in sample #S2 at T=4.2K. 
Inset: the frequency dependence of the real and imaginary parts (open and full circles respectively) 
of the penetration depth Aac(-B,/) from which Ls is deduced. Solid lines are theoretical fit to 
Eq.(l) with Ac = oo, aj^ = lOnOhm.cm and Ls = 14.9/xm. 




FIG. 2. The inverse of phenomenological slippage length ls{B) for a vortex array at a rough 
surface. It is deduced from the Ls data in Fig.l according to the definition Is = Ls/j^qe/B in 
Eq.(l). The effect of vortex interactions in the collective-pinning regime below Bpk is visible as a 
strong reduction in l/ls{B). Solid lines are power- law fits to the data at the transition. 



12 



